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Abstract 

White's density matrix renormalization group (DMRG) method has been ap- 
plied to an S = 1/2 + 1/2 composite-spin model, which can also be considered 
as a two-leg ladder model. By appropriate choices of the coupling constants 
this model allows not only to study how the gap is opened around the gapless 
integrable models, but also to interpolate continuously between models with 
different spin lengths. We have found indications for the existence of several 

different massive phases. 
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I. INTRODUCTION 



The recent discovery of several families of new materials, [0-0] such as (VO)2P2C>7, 
Sr n _ 1 Cu n 2n -i (n = 2, 3, . . .) and La4 +4 „Cu 8+ 2nOi4 +8n , where the spin chains are coupled 
in a special way to form ladders, gave a new impetus to the study of the properties of low 
dimensional magnetic system. This field was already intensively studied both experimentally 
and theoretically in the last decade due to a large extent to the proposal by Haldane 0] that 
antiferromagnetic spin chains described by an isotropic Heisenberg model develop a gap in 
their excitation spectrum for integer spin, while for half-integer spin the spectrum is gapless. 

The spin-ladder models, beside their relevance to these materials, are of special interest 
for theorists PHl5|] because by appropriate choice of the couplings they can describe both 
spin-1/2 and effective spin-1 models, and they are also related to the models proposed to 
understand the behavior of the so-called high-T c materials. Experimental studies |16|JT7|| 
confirm, that two-leg ladders behave like integer-spin models and have finite gap, while 
materials with three-leg ladders have gapless magnetic excitation spectrum. 

A particular way to construct ladder models is to put a composite spin on every site of 



a single chain and to couple the individual spins in various ways. ||lq , |T9|| This model has 
already been studied numerically using the Lanczos algorithm to calculate the low-lying 
energy levels. This method could, however, be applied to relatively short chains only, and 
therefore the conclusions were sometimes contradictory. 



Recently White p(J has proposed a new procedure, the density matrix renormalization 
group (DMRG) method that allows to calculate the energy of low-lying levels and related 
physical quantities on much longer chains. This led to a great progress in the application of 
the finite-size scaling method and allowed to investigate more complex systems in a larger 
parameter space. 

A natural extension of the usual Heisenberg model of spin chains for S > 1/2 is to include 
higher order polynomials of the bilinear exchange term. In the spin-1 case, where biquadratic 
exchange can be taken into account, generically the spectrum remains massive, pTJ but at 
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some special values of the couplings p2|j23[| the gap might disappear and the model can be 
studied by Bethe's ansatz. The opening of the gap around these critical integrable points 
is, however, not quite settled. 

The aim of this paper is to extend the earlier calculations on the composite-spin model 
to longer chains using the DMRG procedure, and to clarify the phase diagram of the corre- 
sponding ladder model. 

The setup of the paper is as follows. In Sec. II we give a short description of the 
composite-spin and ladder models and their relationship to the integrable models for ap- 
propriate choices of the parameters. The DMRG method and the numerical procedures are 
discussed in Sec. III. The results of our numerical calculations are presented in Sec. IV. 
Finally Sec. V contains a brief summary. 

II. THE COMPOSITE-SPIN AND LADDER MODELS 

In a composite-spin model the spin Si at the lattice site i is composed of two or more 
spin operators c^q, (a = 1,2,...). In the most general case the model contains on-site and 
nearest-neighbor interactions among all spin species with different coupling constants. In 
this paper we will focus on a model, where two s = 1/2 spin species, from which an S = 1 
spin can be constructed, are put on every lattice site. 

The model is defined by starting from the most general isotropic spin-1 model, the 
bilinear-biquadratic model, which usually is written in the form 



(1) 



In the composite-spin representation, where Si can be obtained by adding up the two spin- 
1/2 species denoted by &i and fi, Si — a.i + the Hamiltonian takes the form 

H = (cos 9 - - sin 9) (H + + 2 sin 9{H 2 + H 3 ) 
3(N-1) . „ 

+ 4 smg, (2) 
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where N is the number of lattice sites in the chain, and 

= fa " ^+1 + ^ " ^+1] > ( 3a ) 

i 
i 

i 

TC3 = [(^» ' <^+l)(^ ' ^+l) 
i 

+ {pi ■ T i+1 )(fi ■ a i+1 )\ . (3d) 
We will generalize Eq. (Q) and consider the model described by the Hamiltonian 

Ti = Ao Tio + Ai Tii + A2 7^2 + A3 Tiz , (4) 

with arbitrary couplings Aj. 

Alternatively, instead of considering B{ and T{ as spins sitting on the same site, we can 
treat them as sitting on two parallel chains, or on the legs of a ladder, the &i spins on one 
leg and the T{ spins on the other. As shown in Fig. 7i couples spins on the same leg 
only, the others contain inter-leg couplings. Tin, is the usual coupling between spins on the 
same rung, Ti\ couples spins on neighboring rungs of the legs, while TL3 describes four-spin 
couplings on a plaquette. 

Usually the ladder models are constructed to include Ho and TC2 only. When a strong 
ferromagnetic coupling is applied across the rungs (A2 — > —00), the two spins form a triplet 
and the properties of the S = 1 Heisenberg chain are recovered. |||| In another approach 
Barnes et al. allowed for strong antiferromagnetic interchain couplings and treated TLq as 
perturbation. They have shown the existence of a spin gap for any finite interchain coupling. 

In the composite-spin model, on the other hand, Ti\ is also necessarily included. A 
special feature of the model is that for arbitrary values of A2 and A3 the model is invariant 
under the interchange of Ao and Ai, i.e., the energy levels of the full Hamiltonian satisfy 

E(A ,A 1 ,A 2 ,A 3 ) =E(A 1 ,A ,A 2 ,A 3 ). (5) 
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This can be shown by interchanging the o\ and t% spins on every second site. This rela- 
tionship allows to connect the weak- and strong-coupling limits of the model by a duality 
transformation. To show this let us denote by e the energies of the Hamiltonian in which 
the coupling strength of Tio is chosen to be unity. 



E(A ,A 1 ,A 2 ,A 3 ) = Aoe(^,^,^). (6) 

^0 ^0 ^0 



From Eq. (§]) it follows that 



^i A2 A3. ,Aq A2 A3. 



Ao£(y-,f ,^) = A ie (f (7) 
^0 A o A o A l A l A l 

Introducing the couplings Aj = Aj/Ao (i = 1, 2, 3), we get 

e(A 1 ,A 2 ,A 8 ) = A 1 e(i,^ ) ^). (8) 
Ai Ai Ai 

In what follows we will always work with the Hamiltonian in which Ao = 1 and will drop 
the tilde over the couplings. 

For A2 = A3 = this relationship reduces to a usual duality relationship, which connects 
the energies in the < Ai < 1 region to those in 1 < Ai < 00, 

e(A 1 ) = Aie(l/A 1 ). (9) 

When Ai < 0, the above deduced duality relationship connects the lowest-lying levels of 
one region to the highest-lying levels in the other region. Since in the numerical calculations 
a few low-lying levels can only be calculated with sufficient precision, a more useful relation 
can be derived in this case by comparing the energies of Ti defined by Eq. with A = 1 
to that of Ti' = —Ti. Denoting by e'(Xi, A2, A3) the energies of this model, 

e'(\ 1 ,\ 2 ,\ 3 ) = -\ 1 e(±- ,^,^)- (10) 

Ai Ai Ai 

We used these relations to check the accuracy of the numerical calculations. 

According to Eq. (0) the spin-1 model can be reproduced when Tio and Ti\ are taken with 
equal coupling strength, Aq = Ai. When this coupling is antiferromagnetic, it gives rise to 



a valence-bond configuration |24| of the neighboring spins, which is a characteristic feature 
of the Haldane phase. When no direct coupling exists between the &i and t\ spins, nothing 
ensures a priori that these spins appear in the symmetric 5* = 1 configuration only. When, 
however, A2 and A3 are small compared to Ao = Ai, the level structure of the true spin-1 
model and that of the composite-spin model is such that their low-lying parts coincide, thus 
their behavior is similar. Therefore the phase diagram of our ladder model will contain 
phases characteristic to the bilinear-biquadratic model, but also new phases may appear. 

The normalization Ao = 1 is equivalent to considering the region —tt/2 < 9 < tt /2 
only. The bilinear-biquadratic model has a rich phase structure in that region. It has 
been thoroughly investigated both analytically pT| , p5| by using, e.g., the mapping to the 
Wess-Zumino-Witten model, or to the 9-state Potts model and numerically. |26|-p9| The 



region tt/2 < 9 < 3tt/2 is somewhat less interesting, since the model is ferromagnetic for 



tt/2 < 9 < 5tt/A. Beyond that a new phase might appear, as predicted by Chubukov, |31 
although so far its existence has not been confirmed by numerical calculations. |3lJ This 
problem is, however, outside the scope of this paper. 



#tb = — vr/4 corresponds to the Takhtajan-Babujian integrable model with gapless 
excitation spectrum. This point is the critical point of a second-order Ising-type phase 
transition with a gap opening linearly around the transition point. In a chain with periodic 
boundary condition (PBC) the ground state is doubly degenerate for 9 < —it /A producing 
a dimerized phase, while for — tt/A < 9 < tt/4 the Haldane phase appears with a non- 
degenerate singlet ground state. This latter region includes the isotropic Heisenberg point 
at 9 = and the exact nearest neighbor valence-bond state (VBS) |24j at #vbs — arctan(l/3). 



#ls = 7r /4 is another integrable point related to the Lai-Sutherland [^3| model, where 
the gap vanishes again. For 9 > tt/4 a trimerized massless phase appears. p8|-p0| In the 
composite-spin representation these points lie on the line Ai = 1, at A2 = A3 = —4/3 for the 
Takhtajan-Babujian, at A2 = A 3 = for the Heisenberg, at A2 = A3 = 4/5 for the VBS and 
at A 2 = A 3 = 4 for the Lai-Sutherland points. 

Furthermore, the composite-spin Hamiltonian can be transformed into a nonlinear a 
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model and the results obtained in field theory can be related to our model. At Ai = A2 = 
A 3 = 0, when the model reduces to two decoupled spin-1/2 chains, the field theory possesses a 
gapless excitation spectrum because the topological angle is an odd multiple of it. According 
to the field-theoretical calculations |H] the perturbations around this multicritical point are 
generically relevant and give rise to the opening of an energy gap, except for special points 
or lines. The Takhtajan-Babujian and Lai- Sutherland models belong to such exceptions. 
Consequently, it is expected that by varying the Aj parameters the gap vanishes only at the 
critical points and on the phase boundaries between the various phases. 

In order to simplify the calculations, first we study the effect of TCi, Ti.2 and H,3 separately 
for fixed Ao = 1. Then the calculation will be extended to a two-parameter plane, (Ai, A2) 
by choosing A3 = A2. In order to determine the phase diagram we examine the low-lying 
energy levels along different paths connecting those points in the parameter space, where the 
model is integrable. The ground-state configurations will also be investigated by calculating 
the local energy of a bond, the two-point correlation functions and a short-range order 
parameter. 



III. NUMERICAL PROCEDURES 



We have performed numerical calculations by applying the DMRG method p(J on the 
model defined by the Hamiltonian in Eq. (§)-(§)• This is a real-space renormalization method 
where the lattice is built up gradually to the desired length, systematically truncating in 
the meantime the Hilbert space by keeping only the most probable states. 

Since the DMRG method works best for systems with free ends, we will consider our 
composite-spin model with open boundary condition (OBC). An unfavorable consequence 
of OBC for the spin-1 model and consequently for the ladder models as well is that the 
degeneracies in the spectrum may be different from that obtained for a closed ring. In the 
nearest-neighbor valence-bond configuration, e.g., free s — 1/2 spins remain at the ends of 
the chain, giving rise to a fourfold degenerate ground state. In the dimerized state, on the 
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other hand, the twofold degeneracy of the ground state is lifted if the number of sites is an 
even number. This makes the analysis of the spectrum more difficult. 

Moreover, in the case of OBC the total momentum is not a good quantum number. 
Therefore, only the total spin St and its projection to the quantization axis, S T can be used 
to classify the energy levels. Since in the isotropic case the SU (2) symmetry is satisfied, the 
spectrum was analysed by calculating a few low-lying levels of the different S T sectors. St 
was determined from the degeneracy of the levels. 

Because our aim is to determine the overall behavior of the energy spectrum and to 
identify from it the possible phases, in most of the calculations we have used the less accurate 
version of DMRG, the so-called infinite-lattice method. In this algorithm the lattice is 
built up by adding two lattice sites in each step. The results obtained for chains with 
N = 4, 8, 16, 32, 48, 64 sites were used in a finite-size scaling procedure to extrapolate to 
infinite lattice. On the other hand, close to the critical points and in those regions where 
the gap is small, the finite-lattice method was applied with two or three iteration cycles to 
determine the energies more precisely. 

In some cases, before doing the DMRG calculation on long chains, the energy spectrum 
was determined by exact diagonalization on short chains with N = 4,6, 8. This gave an idea 
of the sequence of the levels for different choices of the couplings and allowed to determine 
which states have to be targeted in the DMRG procedure. 

Considering our limited computational resources, we had to restrict rather drastically 
the number of states, M to be kept in the DMRG algorithm. A number, that can be used to 
characterize the numerical accuracy of the DMRG method is the discarded density-matrix 
weight (truncation error). |2(J In the calculation of the ground-state energy for ladders with 
N = 2 x 64 (N = 2 x 32) sites it was largest close to the critical points with values of the 
order 1CT 5 — 1CT 6 (1CT 7 — 1CT 8 ), while around the VBS point, where the gap is large, it was 
as small as 1CT 10 — 1CT 11 (1CT 11 — 1CT 12 ). For the excited states the truncation error was 
worse by one order of magnitude. 

The real error in the DMRG procedure can, however, be much worse. |33|] It was estimated 
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by comparing the energies obtained for the ladder model using the DMRG procedure to those 
of the spin-1 bilinear-biquadratic model or the spin-1/2 chain, where much better accuracy 
can be achieved. For chains with N = 2 x 32 sites an agreement up to three or four decimal 
places has been achieved after the second iteration cycle of the finite-lattice method, if M 
was chosen to be M = 64, while for chains with N = 2 x 64 sites the same accuracy required 
to retain M = 76 or 84 states. At the VBS point the infinite-lattice algorithm without the 
iteration cycles provided almost the same accuracy. 

In the extrapolation procedure the N — > oo limit of the gap was determined by fitting a 
form AE(N) = A + a/N or AE(N) = A + b/N 2 to the energy differences measured between 
the energy levels of the various St sectors. In principle the first has to be used when the 
gap vanishes, A = 0, while the second gives the correct asymptotic behavior when the gap 
is finite. |34j] When the gap is finite but the chain length is not long enough to observe the 
correct parabolic behavior, the linear fit gives a lower-bound estimate. |35[ We have used 
the parabolic fit, whenever deviation from the linear dependence on 1/N was observed. 



IV. NUMERICAL RESULTS 

In this section we present our numerical results for the low-lying levels of the model for 
various choices of the couplings. For a finite ladder the ground state is always a spin singlet. 
Above this level there are two low-lying triplet excitations, which at Ai = A 2 = A3 = 
have equal energy. At this point, they correspond to independent excitations on the two 
decoupled legs. These levels are followed by several 5V = 0, 1, 2, ... , excitations. The energy 
differences between all these levels scale as 1/N, indicating a gapless excitation spectrum. 

For finite values of Ai, A2 and A3 the degeneracy of the two low-lying triplet levels is in 
general lifted. This is shown along the Ai axis in Fig. 0. In a finite chain with N sites 
the energy difference between the lowest triplet and the ground-state singlet levels will be 
denoted by AEiq(N). When this quantity scales to a finite value, the ground state is a non- 
degenerate singlet with a finite gap. When, however, it scales to zero, we have to consider 
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the next level. This is either the second triplet, or it may have been crossed by the lowest 
St = 2 level. The energy difference between this quintuplet and the ground-state singlet 
levels will be denoted by AE2o{N), while AE2i(N) denotes the energy difference between 
this quintuplet and the lowest triplet levels. 



A. Gap opening due to Hi 

Let us consider first the effect of TCx, since earlier calculations by exact diagonaliza- 
tion on relatively short chains could not determine satisfactorily the way the gap is opened. 
For Ai > the energy of the lowest triplet states and that of the singlet ground state 
come exponentially close to each other, and become degenerate in the thermodynamic limit 
(N — » oo). This gives the well known fourfold degeneracy of the VBS-like state in a finite 
chain with OPB. 

Therefore the relevant energy gap can be obtained most easily by measuring the energy 
difference AE2i(N) between the lowest-lying levels of the St = 2 and St — 1 sectors. This 
energy difference as a function of Ai is shown in Fig. [5] for Ai > 0. The lines connecting 
the values calculated at several distinct points are only guides to the eye. The inset in the 
figure shows in more detail the behavior of the gap close to the integrable point. 

In the extrapolation to N — > oo a 1/N fit was applied close to the critical point, because 
the chain lengths were not long enough to observe the parabolic behavior. As an indication 
of the accuracy of our calculation we mention, that at Ai = 1, where the Haldane gap should 
be recovered, we obtained A21 = 0.41(1), in reasonable agreement to three decimal places 
with its best estimate. [f34| , |36f The ground-state energy divided by the chain length converges 



to Eq/N = 1.401484(0), in agreement to six decimal places with earlier calculations. [36| 

As a further check we have also performed calculations for large values of Ai and used 
the self-duality relationship in Eq. (|9|) to obtain the gap for weak couplings. The results 
obtained in this way are shown in the inset. From this we could also conclude that our 
finite-size calculations are correct to at least 4 decimal places. Even though the accuracy of 
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our calculations is limited, we have found a linearly opening energy gap for Ai > 0. This is 
the first calculation where this could be demonstrated numerically 

For negative values of Ai the coupling between the chains is ferromagnetic, while the 'en- 
chain' coupling remains antiferromagnetic. The low-lying part of the spectrum is essentially 
different from that of the Ai > regime. The degeneracy of the two St — 1 triplet levels 
is again lifted, but even the lowest of them will remain separated from the singlet ground 
state. So the energy difference to be studied is AE W (N). This quantity is shown on the 
left hand side of Fig. || for Ai < 0. As it can be seen, the ground state is a non-degenerate 
singlet with a small, but finite gap to the lowest lying excitations. 

Since the gap is small, the value of A 10 has been calculated using the finite-lattice method. 
This procedure with three iteration cycles and a 1/N fit gave at Ai = —1, e.g., A 10 = 0.11(4) 
as a lower-bound estimate for the singlet-triplet gap. Furthermore, the lowest quintuplet 
level has been found to become degenerate with the lowest triplet state in the thermodynamic 
limit. Thus a gap separates the ground state from the continuum of excitations. 

For small values of Ai the accuracy of the calculations has been checked using the rela- 
tionship in Eq. ([TD|) . The gap is found to open linearly. Thus Hi is relevant for both signs 
of A x . 

On the other hand, for large negative values, in the region Ai < —10, we have found the 
vanishing of both A 10 and A 20 . This indicates that a new phase with gapless spectrum may 
appear there. However, due to the smallness of the gap around Ai = —1 we were unable to 
locate the point where the transition occurs. 

B. The effect of H 2 and H 3 

The effect of the inter-leg coupling "Hi on the decoupled chains described by 7io for 
Ai = A3 = has been considered by several authors. [|5||8|,|T0"1| Hida found that the gap is 
finite for A2 < A 2c = —0.6, while Dagotto et al. showed some evidence that the critical 
coupling is closer to zero, perhaps A 2c ~ —0.4. On the other hand Watanabe et al. || and 
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Hsu and Angles d'Auriac |10| argued that the interchain coupling defined by TC2 is always 
relevant and A 2c = 0. 

We repeated some of these calculations to confirm that this coupling is in fact relevant 
for both signs of the coupling A 2 . For A 2 > the ground state is a non-degenerate singlet 
even in the thermodynamic limit with a rather large gap. At A 2 = 4/3, e.g., a singlet-triplet 
gap with A10 = 0.72(9), and a singlet-quintuplet gap with A 20 = 1.43(4) has been found. 

For A 2 < the spectrum is similar to that for Ai > 0, where the ground state becomes 
fourfold degenerate in the thermodynamic limit. The gap, however, is rather small. At 
A 2 = —4/3 the calculation with the finite-lattice method and a 1/N fit gave A 2 i = 0.11(3), 
in agreement with previous results. || 

The plaquette coupling TC3 also turned out to be relevant for both signs of the coupling, 
but in some sense its effect is opposite to that of H.2- A fourfold degenerate ground state is 
obtained with a finite quintuplet-singlet gap for A3 > 0. At A3 = 4/3, e.g., A 20 = 0.28(5) has 
been found. For A 3 < the ground state remains singlet. At A 3 = —4/3 the triplet-singlet 
gap is A 10 = 0.26(8). 

Because of this opposing effect of these couplings, next we considered the competition of 
H.2 and Ti^ by choosing A 2 = A3. For negative values of A 2 the low-lying part of the energy 
spectrum resembles very much that found above for small negative Ai. A finite gap develops 
between the singlet and triplet levels as shown in Fig. [|. At A 2 = —4/3 the finite-lattice 
algorithm and a 1/N fit gave A 10 = 0.17(5). Close to the critical point the extrapolated 
lower-bound values of the gap are even smaller, therefore the available chain lengths are still 
not long enough to determine the character of the opening of the gap. 

On the other hand, for A 2 = A3 > the analysis of the low-lying energy spectrum has 
shown that both A 10 and A 20 scale to zero, i.e. the spectrum is gapless. Since the system 
remains critical in an extended region for A 2 = A3 > 0, a Kosterlitz-Thouless-like transition 
may occur at Ai = A 2 = A3 = to the massive phase at A 2 = A3 < 0. This could explain 
the slow opening of the gap in that region. 

These results are only partially consistent with the field theoretical predictions, which 
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states that the spin-1/2 integrable point is unstable against all perturbations. Our calcu- 
lations indicate that although both 7i 2 and are relevant operators in the ladder model 
for both signs of the couplings A2 and A3, the model remains critical at least along the 
A 2 = A 3 > half line. 

To summarize our findings, three types of the spectrum have been found. An asymp- 
totically fourfold degenerate ground state, charateristic of a VBS-like state is obtained for 
Ai > 0, A2 = A3 = 0, for A2 < 0, Ai = A3 = and for A3 > 0, Ai = A2 = 0. A truely non- 
degenerate singlet ground state is found for Ai < 0, A2 = A3 = 0, for A2 > 0, Ai = A3 = 0, 
for A3 < 0, Ai = A2 = and also for A2 = A3 < 0, Ai = 0. Finally the spectrum is gapless 
for A 2 = A 3 > 0, Ai = 0. 

C. Competition of Hi, TL2 and H3 

After having determined the spectrum along the Aj axes and the A2 = A3 line, next 
we consider the phase diagram in the parameter space spanned by Aj. To simplify the 
calculations we restrict ourselves in the remaining part of the paper to a two-parameter 
plane by choosing A 2 = A 3 . We will study especially the neighborhood of the Ai = 1 line, 
which corresponds to the bilinear-biquadratic Hamiltonian. 

In the Takhtajan-Babujian point of the bilinear-biquadratic model, which corresponds to 
Ai = 1, A2 = A3 = —4/3 in our model, the spectrum is gapless in the thermodynamic limit. 
Both the singlet-triplet (AE W (N)) and the triplet-quintuplet (AE 2 i(N)) energy differences 
vanish as 1/N. They behave, however, quite differently as we move away from the critical 
point. This is shown in Figs. |5| and || 

For A2 > —4/3 the lowest triplet level becomes asymptotically degenerate with the 
ground-state singlet, so the relevant gap is between the triplet and quintuplet levels (A21). 
This level structure is the same as along the Ai > line, thus a VBS-like state is obtained 
in this part of the phase space. The opening of the gap is very slow, even our longest chains 
are too short to obtain a reliable estimate of the gap close to the critical point. The error 
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of the extrapolated value close to the critical point is indicated in Fig. || by the size of the 
symbol. 

For A2 < —4/3 a finite gap develops between the ground-state singlet and the lowes 
triplet levels, in the same way as for A2 = A3 < 0. A lower-bound estimate of the gap 
obtained by a 1/N fit is shown in Fig. |^. The results are in agreement with a linearly 
opening extrapolated gap Ai . 

It is known, however, that the A2 < —4/3 region corresponds to a dimerized phase, where 
instead of a singlet ground state a doubly degenerate ground state would have to be found. 
Therefore, we have considered several higher lying levels to search for the other singlet level 
that would become degenerate with the ground state. We have not found any such level. 
This can be understood by recalling that in the case of OPB the dimerized phase gives a 
truely twofold degenerate ground state only if the number of sites is odd, i.e., the number 
of bonds is even. For chains with even number of sites the energy of the dimerized state 
depends on whether the bonds at the ends are strong or weak. Inspection of the low-lying 
levels is therefore not sufficient to distinguish a real non-degenerate singlet ground state from 
a dimerized state. We will return to this problem later, when the dimer order parameter 
will be discussed. 

One of the two kinds of behavior found for A2 = A3 > — 4/3 and A2 = A3 < —4/3 appears 
whenever one moves away from the Takhtajan-Babujian integrable point in any direction 
in the (Ai,A 2 = A 3 ) plane. Varying, e.g., Ai around Ai = 1 at A 2 = A 3 = —4/3, a non- 
degenerate singlet ground state is found for Ai < 1, while for Ai > 1 the singlet-triplet gap 
disappears forming a fourfold degenerate ground state. This result is in agreement with the 
assumption that the Takhtajan-Babujian point is generically unstable against perturbations, 
except along the phase boundaries. 

Our calculation is, however, not accurate enough to locate this critical line in the phase 
space that separates the two kinds of behavior. In the schematic phase diagram shown in Fig. 
[7|the boundaries are therefore indicated by wavy lines. That the phase boundary between the 
VBS-like Haldane phase and the supposedly dimerized phase connects the Ai = A2 = A3 = 
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and the Takhtajan-Babujian points has been confirmed by looking at the energy spectrum 
of our model along the trajectory parametrized by Ai = 1 — A, A 2 = A 3 = — 4A/3, where A 
varies between and 1. For small values of A a fourfold degenerate ground state is recovered, 
while for A close to unity a non-degenerate singlet ground state is found. 

As mentioned before, this singlet ground state was found along the A x < 0, A 2 = A 3 = 
line as well. Whether this axis also belong to the dimerized phase or not, cannot be 
determined from the low-lying spectrum alone. We will return to this problem later. 

Extending the calculations to A 2 = A3 > on the line Ai = 1 the VBS-like state survives 
in a finite range, including the point A 2 = A3 = 4/5, where the exact nearest-neighbor valence 
bond state is recovered. At this point the lowest singlet and triplet levels are degenerate for 
any finite chain lengths. The gap reaches its maximum value with A 21 = 0.8404(7), which 
after including the appropriate scaling factors due to our normalization to Ao = 1, agrees to 
5 digits with the known result. [|37j1 

For larger values of A 2 along the Ai = 1 line the Lai- Sutherland model and the trimerized 
phase of the bilinear-biquadratic model cannot be reached by this ladder model. The extra 
levels introduced by the composite-spin representation will be low lying and will lead to a 
gapless new phase, as has already been found for Ai = 0. This massless phase has been 
found for large positive A 2 = A3 values also for Ai < 0. Thus it is stable in an extended 
range of the couplings. 

D. Dimer order 

Even though the structure of the low-lying part of the energy spectrum indicates the 
existence of various phases, its knowledge, as discussed above, is not sufficient to clarify 
unambiguously the character of the ground state. We have, therefore, calculated several 
quantities in the ground state, like the local magnetization (Si) = (ff{) + (7^), the two-point 
correlation function (SiSj) = ((ai + ri) (aj + fj)) and the local energy E\ oc = (H(i, i+ 1)), 
where Ti(i, i + 1) contains the couplings between spins on sites i and i + 1. 
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The correlation function falls off exponentially both in the Haldane and the dimerized 
phases, turning to power-law like on the phase boundary, but the chains are still too short 
to distinguish clearly between these two possibilities. A better procedure could be to look 
at the local magnetization in the lowest triplet state. Due to the free end spins of the exact 
nearest-neighbor valence-bond configuration the local magnetization is finite close to the 
chain ends in the VBS-like Haldane phase. Approaching the phase boundary, the extra 
spin becomes less and less localized and it spreads out homogeneously in the dimerized and 
gapless phases. 

The most useful procedure is, however, to consider the so-called short-range dimer order 
parameter. It can be defined by taking the difference of the local energy on neighboring 
bonds in the middle of the ladder, 

77 (H(i,i + l))-(H(i+l,i + 2)) .... 

Ksido \ [{H(i, i + 1)> + (H(i + + 2)>] • [LL) 

It is expected to have different behavior if the ground state is unique, fourfold degenerate, 
or twofold degenerate as in the dimerized state with spontaneously broken translational 
symmetry. 

A simpler quantity can be defined by taking the bilinear part of the coupling only 

(Si ■ Si+i) — (Si + \ ■ Si+2) 



S 



srdo 



(12) 



(Si ■ Si + \) + (Si+i ■ Si+2) 
In Fig. § we present our results for the short-range dimer order parameter S^do measured 
in the middle of the ladder as a function of the inverse of the chain length for a few points of 
the (Ai, A2 = A3) phase space. The parameter -R sr do not shown in the figure gives the same 
kind of behavior. 

There is clearly an extended region in the parameter space, where the short-range order 
parameter scales to zero in the thermodynamic limit. At the VBS point itself SWdo was 
found to be zero for short finite chains already. The vanishing of S sr do happens not only in 
the Haldane phase, but also along the negative Ai axis and in a neighborhood of it, where 
the spectrum was undistinguishable from the spectrum of the supposedly dimerized state. 
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For large negative values of A2, however, both for positive and negative values of Ai, a 
transition was found to the state where -R S rdo and S ST ^ have non-zero value and the local 
magnetization at the ends of the chain vanishes. Thus in the region where the ground state 
was found to be a non-degenerate singlet, the short-range dimer order allows to distinguish 
two regimes. In the region denoted as 'dimerized phase' in Fig. |7|, both -R sr do and S ST do scale 
to a finite value. In this phase the local magnetization at the ends of the chain is also absent. 
In the other region marked as 'massive', the dimer order disappears. The boundary between 
them is drawn schematically only in Fig. [7|, its exact shape could not be determined, except 
that it has to go through the Ai = A2 = A3 = point. 

V. CONCLUSIONS 

In the present paper we have considered a two-leg ladder model constructed from a 
composite-spin model. Beside the usual Heisenberg coupling between the spins on the legs 
inter-leg coupling between spins on the same and neighboring rungs have been introduced, 
as well as four-spin plaquette couplings. 

For special values of these couplings this ladder model is equivalent to the spin-1/2 
Heisenberg model or the spin-1 bilinear-biquadratic model in the sense that the low-lying 
parts of the spectra are identical. Thus massless, VBS-like and dimerized phases are expected 
to appear, but the phase diagram can be even richer. 

The behavior of the energy differences between low-lying levels has been calculated using 
the DMRG method and the gap extrapolated to the thermodynamic limit (N — > 00) was 
obtained with the method of finite-size scaling. We have also considered the short-range 
dimer order parameter. 

It has been found that the spin-1/2 Heisenberg ladder is unstable against most of the 
perturbations that couple the two legs, but there is a small range of the parameters where 
the combination of the perturbing operators is irrelevant. 

Four different kinds of behavior was observed, as shown in Fig. |7|. In a region along and 
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near to the line Ai = 0, A2 = A3 > the spectrum remains gapless. Everywhere else a gap 
is developed in the spectrum. Three regions can, however, be distinguished depending on 
whether the ground state is non-degenerate, twofold or fourfold degenerate. The last case 
is easily detected by studying the asymptotic degeneracy of the spectrum. This phase is 
present around the Ai = 1, A 2 = A 3 = 4/5 point, where the exact nearest-neighbor valence- 
bond state is recovered. The extra s = 1/2 spins at the ends of the chains become more 
and more delocalized as we move away from the VBS point and disappear at the phase 
boundary, where the transition is either to the gapless phase or to a dimerized one. 

The Takhtajan-Babujian point at Ai = 1, A2 = A3 = —4/3 lies on the phase boundary 
between the VBS-like and the dimerized phases. The transition here could be detected not 
only in the change of the character of the spectrum, but also by the appearance of a finite 
short-range dimer order parameter. This quantity was used to distinguish the dimer phase 
also from the 'massive' phase with non-degenerate singlet ground state, where the dimer 
order disappears again, in fact much faster than 1/N. 

Unfortunately, our limited computational resources did not allow us to determine the 
precise location of the critical lines separating the various phases. 
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FIGURES 

FIG. 1. Schematic plot of the spin couplings in Ho, Hi, H2 and H3 between the spins c?j and 
Ti on the two legs of a ladder. 

FIG. 2. Low-lying energy spectrum of the Hamiltonian as a function of Ai (A2 = A3 = 0) for 
chain length N = 6. 

FIG. 3. The energy difference between the low-lying levels of our model as a function of Ai at 
A 2 = A 3 = for chain lengths 8 < N < 64. For Ai > A£ 2 i, while for Ai < AE W is shown. The 
dotted line is the gap extrapolated for N — ► 00. The inset shows the behavior for small values of 
Al 

FIG. 4. The energy difference AEiq as a function of A2 at Ai = A3 = . The symbol x shows 
the extrapolated value of the gap for A 2 = —4/3. 

FIG. 5. The quintuplet-triplet energy difference AE2i{N) as a function of A 2 = A3 at Ai = 1. 
Close to the critical point the size of the symbol for the extrapolated gap indicates the error. 

FIG. 6. The same as Fig. 5 for the triplet-singlet energy difference AEio(N). 

FIG. 7. Schematic phase diagram of the model in the (Ai,A 2 = A3) plane. The wavy lines 
between the four different phases indicates that the location of the phase boundaries has not been 
determined accurately. 

FIG. 8. The short-range dimer order parameter S sr d calculated in the middle of the chain 
versus the inverse of the chain length at a few points of the phase space (Ai, A 2 = A3). 
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